The analytical inr scattering amplitude and the light scalars 

N.N. Achasov a and A.V. Kiselev a ' 6 
a Laboratory of Theoretical Physics, Sobolev Institute 
for Mathematics, 630090, Novosibirsk, Russia 
b Novosibirsk State University, 630090, Novosibirsk, Russia 
(Dated: March 17, 2011) 

Abstract 

In this work we construct the -kit scattering amplitude Tq with regular analytical properties in 
the s complex plane, which describes simultaneously the data on the tttt scattering, <j) ~^ 
decay, and tttt — > KK reaction. The chiral shielding of the <r(600) meson and its mixing with the 
/o(980) meson are also taken into account. The data agrees with the four-quark nature of the 
cr(600) and /o(980) mesons. The amplitude in the range — 5m^ < s < 0.64 GeV 2 also agrees with 
results, obtained on the base of the chiral expansion, dispersion relations, and the Roy equations. 

PACS numbers: 12.39.-x 13.40.Hq 13.66.Bc 
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I. INTRODUCTION 



Study of light scalar resonances is one of the central problems of nonperturbative QCD, 
it is important for understanding both the confinement physics and the chiral symmetry 
realization way in the low energy region. The commonly suggested nonet of light scalar 



mesons is / Q (600) [or cx(600)], K*(800) [or k(800)], /o(980), and a (980) Light scalar 
mesons are intensively studied theoretically and experimentally in different reactions. 

In Refs. |2| we described the high-statistical KLOE data on the <fi — > 7r°7r°7 decay 
[3I simultaneously with the data on the rrn scattering and the ttty — > KK reaction. The 
description was carried out taking into account the chiral shielding of the <r(600) meson (4, ^ 
and its mixing with the /o(980) meson. It was shown that the data do not contradict the 
existence of the cr(600) meson and yield evidence in favor of the four-quark nature of the 
<r(600) and / (980) mesons. 

This description revealed new goals. The point is that at the same time it was calculated 
in Ref. js] the wir scattering amplitude in the s complex plane, basing on chiral expansion, 
dispersion relations, and Roy equations. In particular, the pole was obtained at s — M% = 
(6.2 — 12.3i) m^, where 

M a = AAl+l 6 - z272i? 2 5 MeV , (1) 

which was assigned to the a resonance. 

Aiming the comparison of the results of Refs. Q and it is necessary to build the nn 
scattering amplitude with correct analytical properties in the complex s plane. The point 
is that in Ref. (2] S matrix of the tctt scattering is the product of the "resonance" and 
"background" parts: 



S-mr 'S'feacfc S res , (2) 

and the S res had correct analytical properties, while analytical properties of the Sback in the 
whole complex s plane were not essential for the aims of [2|, where the physical region was 
investigated, and Adler zero existence |?J together with the poles absence on the real axis of 
the s complex plane were demanded. 

In this paper we present the tttt scattering amplitude with correct analytical properties 



in the complex s plane and the data description obtained with this amplitude 



The 



comparison with the results of Ref. [6| is also presented. 

All formulas for the <fi — > (5*7 + p°7r°) — > 7r°7r°7 reaction [S = /o(980) + a(600)] are shown 
in Sec. [TT1 Our new parametrization of the background amplitude is presented in Sees. Ill II 
and IIVI The results of the data analysis are presented in Sec. |V] A brief summary is given 
in Sec. ED 



II. THEORETICAL DESCRIPTION OF THE <j> -»• (/ (980) + <r(600))7 -> jtt°tt AND 



(/) -»■ A -> 7^° REACTIONS 

In Refs. (9), [lo| it was shown that the dominant background process is </> — > Tt°p — > 77r°7r°, 
while the reactions e + e~ — > p — > tt°u — > 77r°7r° and e + e~ — > to — > n°p — > 77r°7r° have a 
small effect on e + e~ —>■</>—>■ 77r°7r° in the region m„o n o = m > 900 MeV. In Ref. [ll| it 
was shown that the <j) ~~** 7r °P ~~ ► 77r°7r° background is small in comparison with the signal 
-> 7/o (980) ->■ 77T 7r° at m > 700 MeV. 

The amplitude of the background decay 4>{p) — > n°p — > 7(g)7r (/ci)7r°(/c 2 ) has the following 
form: 



Mb ac k — Fbe % g p7r o^ ) g p7r o 7 (j) a p u €sqe€a^fj,i>^Su 



kl/i^uj 



+ 



ko„.k 



(3) 



.D p (q + k 2 ) DJq + h 
Here, constants Ff, and 5 take into account pn rescattering effects [12]. Note that in this 



work and our previous work it was assumed that Ff, = 
In the K + K~ loop model, -> K + K~ -> 7(/ + cr) 
amplitude of the signal — >■ 7(/ + cr) — > 77r°7r° is 



a 



9h11|. above the KK threshold the 



M, 



sty 



g{m) 



W-aM)T(rr^V)xi6,, 

[jpq) J V / 



(4) 



where the K + K — > n tt amplitude, taking into account the mixing of /o and a mesons, 



T (K + K- 



7rV 



E 

R,R' 



9RK+K- GRR'9R'n°Tr° 



167T 



where R, R' = fo, a, 



Ob = o R +0 



KK 
B ~<~ U B ) 



(5) 



(6) 



where c^ 71 " and 5^ are phases of the elastic background of the im and KK scattering, 



respectively, see Refs. 



14 



rzj. 
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Note that the additional phase 5^ changes the modulus of the KK — > n n amplitude 
under the KK threshold, at m < 2m K . Let us define 



5 kk 



K 



I analytical continuation of e 



i6§ K 



m > 2i7iK ; 
m < 2rriK ■ 



Note also that the phase 5 B n was defined as 5b in Refs. 
The matrix of the inverse propagators [10] is 

/ 

Grr' = G RR > (m) = 



id, 



ill. 



V 



D f ( m ) -^fA m ) 

•n /o(T (m) D a (m) 



(7) 



n 



f0T 



m 



E 



0W6 -n^(m) + c foa , 



a,b 9f ab 

where the constant C/ 0(J incorporates the subtraction constant for the transition /o(980) — > 
(0~0~) — > (j(600) and effectively takes into account the contribution of multiparticle inter- 



mediate states to fo 4-> a transition, see Ref. [10j. The inverse propagator of the R scalar 



meson is also presented in Refs. 



14 



23j: 



D R {m) = m 2 R - m 2 + J2[ReU a R b (m R ) - Uf(m% (8) 

ab 

where ^ a! \ReTl R h (m R ) — U R b (m 2 )] = ReU R (m R ) — U R (m 2 ) takes into account the finite width 
corrections of the resonance which are the one loop contribution to the self-energy of the R 
resonance from the two-particle intermediate ab states. 

For pseudoscalar a, b mesons and m a > nit,, m > m + one has 



m_ <m< m + 



9Rab 
167T 



m + m_ , mi, 
9 In— + 



1 \jm* — m"_ — \/m* — m 
+p ab [i + -In _ 

" Jm z — m_ + Jm 



ml 



(9) 



9Rab 
167T 



77- m — 



m a 



- \Pab(m)\ + 



+ — \pab{ m )\ arctan 

7T 



m 4 



in 2 — m 2 _ 



(10) 
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m < m_ 



9Rab 

16n 



— ln 

Trm z m„ 



—Pab(™>) ln 
7T 



mi — m? — \ ml — m 2 



ml — m 2 + */ ml — m 2 



AM™) = 

The constants are related to the width 



9 2 

mi ml 

1 -)(1 -) , m + = m a ±m b 

m z wr 



T R (m) =Y,r(R ^ ab,m) =J2 



9Rab 

167rm 



Pab(m) 



(11) 



(12) 



(13) 



ab ab 

Note that we take into account intermediate states irir, KK,r)r),r)'r),r)'r)' in the /o(980) 
and a (600) propagators: 



% = n? + - + n£*° + nf + uf^° + nf o + + (14) 

and also for the a(600). We use g^ K o^o = gf K+K-, 9^°^° — fl , / 7r+7r-/v / 2, the same for the 
a(600), too. 

For other coupling constants the naive four-quark model predicts 



9 form - -9fov'v' 



2^2 



V2 



3 9f K+K-, 9fov'V — £-9foK+K-', 



9or\r] 9ur]r]' 



V2 



?<T7r+7r _ ) 9or]'r\ 



3y/2' 



The definition of fi^o^o, gn m , gR V >rf takes into account the identity of the particles. As 
these relations are approximate, we introduce the effective correction coefficients x a and Xf : 



9 farm ~ 9fov'v' 



2^2 



9f K+K- Xfo, 9 for)' r, 



V2 



■ 9f K+K- Xf ; 



V2 



gcrriri garm' g g<nr+ir ^aj g<Jr]'r}' 3y^2 ^ <T7T+n ~'^ a 



In the K + K loop model g(m) has the following forms (see Refs 



22 



23|). 
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For m < 2m K + 



g(m) 



2(2tt) 



e f 1 — p 2 (m 2 ) 

HoK K \ 1 + ~ ~ 7T7~ ^r - x 



p 2 (mV) — p 2 (m 2 ) 



2\p(m 2 )\ arctan^— — - p(m 2 )A(mJ) + mp{m\) - 



/ 1 

arctan 

V b( m ) 



2\ n 



where 



T- 



4ml 



/>('» )-- v • 
V 



A(m 2 ) = In 



1 + p(m 2 ) e 2 
1 — p{m 2 ) ' 47r 



For m > 2m K + 



g(m) 



-,9<t>K+K- { 1 



1 - p (m ) 



2(2tt) 2 ^'" f ' p 2 (m 2 ) -p 2 (m 2 ) 



X 



x p(m )(A(m ) - i7r) - p(m^)(A(m^) - ztt) - 
1(1 - p 2 (mj)) ((tt + a(mj)) 2 - (tt + *A(m 2 )) 2 ) 
The mass spectrum of the reaction is 



r(0 ->• 7r°7r°7) _ dT s dT hack {m) dT int (m) 



dm dm dm 

where the signal contribution — > Sj — > 7r°7r°7 



dm 



dT s \PK\\g(m)\ 2 Jm 2 -Aml(ml-m 2 ) 



^2 9RK+K-G RR/gR'TT 
R,R> 



dm 3(47r) 3 m^, 
The mass spectrum of the background process — > pn° — > 7r°7r°7 



' "~~2 256^1 J^dxA^^x), 



dm 



where 



A back {m, x) = ^Y1 \ M back\ 2 
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r <7^ 7 j (ml + 2m 2 m 4 m 2 - 4m^m^ + 2m 4 m 4 p - 



24 ^<PP?r- 

— 4m 2 m 2 m 4 , + 6m 4 m 4 + 2m 2 m % p — 4m 2 m 6 p + m p — 2m^.m 2 — 
-2m 2 mlm 2 m 2 + 2m 4 m 2 p m 2 t> - 2m 2 m A p m\ + 2m 2 m A p m 2 (j) - 2m 6 p m 2 + 

-2ml + 2mJ - m 2 ) {2m 2 ml + 2m> 2 - m 4 ) 1 + 

\ LJ (A m P )\ 

+2Ref - — 1 ) (ml - m 6 m 2 + 2m 4 m 2 m 2 + 

+2m 2 m 1T m 2 — 4m®m 2 — 4m 2 m 2 m 4 + 6m 4 m p + 

+2m 2 m 6 p — 4m 2 m p + m 8 p + m 2 m\m 2 s — 2m\m 2 s + 2m 4 m 2 m 2 b — 

A 2 2-2 2 i o 4-2 2 2-4 2 , r> 2-4 2 r> ~ 6 2 

-4m m 7r m p m^ ) + 2m 7T m p m (j> - m m p m^ + 2m 7T m p m (j> - 2m p m (j> - 

4 4 2-2 4,o 2-2 4,-4 L 

—m^m^ — m m p m^ + 2m 7T m p m ( p + m p m ( 



- 2 2 i - m 2 ) r 4^2" 

m p * 2 = m 2 + 2m 2 — m 2 — m 2 . 



The interference between signal and background processes accounts for 



dm ~ V2 2^m% L^Ant^x) , 



where 



A int {m, x) = -{ml - m 2 )Re J2 M f Ml ack 



167T 



F b Re\e lS g{m) g<Pp7r g p7T1 T° (k + K~ 7r°7r°) 



(m 2 — m 2 ) 2 m 2 . — (m 2 — m 2 ) 2 m 2 
D*{m p ) 



+ 



(m* 2 — m 2 ) 2 m 2 — (m 2 — m 2 ) 2 m* 2 
D* p {m* p ) 



yRejp^e^V 5 #(m)^^o T ^^ 9rk+k-G RR ,g Rl ^^ x 



x 



2\2™*2 



r (m 2 — m 2 ) 2 ??^ — (m| — m 2 ) 2 m 2 {m* 2 — m 2 ) 2 m 2 — (m 2 — m 2 ) 2 m p 



D*{m p ) 



+ 



D*{m* p ) 
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The factor 1/2 in Eq. ( 12"U|) and the factor l/y/2 in Eq. f[2"3"j) take into account the identity 
of pions. 



The S-wave amplitude Tq of the im scattering with 1=0 [lOj, ll5Ml7j is 



rpV, 

1 n 



v y< - 1 



2i8' n n i i 



2ip 7r7r (m) 2ip mr (m 



16tt 



Here 77° = f?o( m ) i s the inelasticity, 77$ = 1 for m < 2m K +, and 



(25) 



<5 u = 5 u (m) = <5™(m) + 5 res (m) 



(26) 



where 5g 7r = d B *(m) (5b in Ref. [lOj) is the phase of the elastic background [see Eq. [6], and 
8 res {m) is the resonance scattering phase, 



M 167F 



gr fl7r7r = ^3/2 gRTT+Tr-. The chiral shielding phase 8^(171), motivated by the a mode 
and desired analytical properties, is taken in more complicated form than in Ref. 
Sec. ED 

The phase 8g = S§ (m) is parametrized in the following way: 



(27) 



sec 



tantff K = f K (m 2 )^m 2 -Am 2 K+ = 2p K f K (m 2 



(28) 



and 



M*** 1 +i2p K f K (m 2 ) 
e B 



1 - i2p K f K (m 2 ) 



(29) 



Actually, e 2 * 5 5' r ( m ) has a pole at 



m = m 



< Trig < 4m 2 , which is compensated by 



the zero in e 2t5 B K ( m ) to ensure a regular ii'i^ — > tttt amplitude and, consequently, the — > 
K + K~ — > 7T7T7 amplitude at < m 2 < 4m 2 . This requirement leads to 

1 1 



fK(m ) 



As in Refs. [2J, for /^(m 2 ) we used the form 



4m|- + - m 2 , 2m K + 



(30) 



, /in. — //ti 

arctanf 3-^ 



A 



K 



(31) 



The inverse propagator of the p meson has the following expression: 



D p {m) 



m, 



2 ■ 2 9pmr 

m — im ' 



48tt 



1 - 



4m 2 \ 3/2 



(32) 



The coupling constants g^K+K- = 4.376 ± 0.074 and g^ p7T = 0.814 ± 0.018 GeV -1 are 
taken from the most precise measurement J24j . To obtain the coupling constant g p7r o 7 we 
used the data of the experiments 25| and [26j on the p — > 7r°7 decay and the expression 



T{p -> 7T 7 ) 



n 2 



P"" u 7 I 2 2\3 



967rm 3 , 



(33) 



the result g p7r o 7 = 0.26 ± 0.02 GeV 1 is the weighed average of these experiments. 



III. THE BACKGROUND PHASE 5% 



The proper analytical properties of the tttt scattering amplitude are two cuts in the s- 
complex plane, Adler zero in Tq [2^] , absence of poles on the physical sheet of the Riemannian 
surface, <r(600) and /o(980) poles in the resonance amplitude on the second sheet of the 
Riemannian surface, and absence of poles on the second sheet in the background amplitude 
in the region 4m 2 < Re(s) < (1.2 GeV) 2 . This applies curtain restrictions on the 5]^. 

Let us represent S 7 ^ in the physical region s = m? > 4m 2 as 



and 



where 



tan(^) 



lm(P nl (s)P n2 {s)) 
Re(P 7rl (s)P 7r2 (s)) 



OjSjTTir 

e B 



C^back cutback 



PZMP*As) Pm(s - ie)P^{s - it) 



i^i(s)P ff2 (s) P 7rl (s + ie)P 7T2 (s + ie) 



(34) 
(35) 



-Ptti(s) = ax - a 2 



Ami 



U n7r (s) + a 3 n 7r7r (4m 2 - s) - a 4 Qi(s) , 



(36) 



71 J Ami S' 



Ami s' — s — is 



(37) 



Kils) 



LAs) 



D^Aml - s)D 2 (Aml - s)D 3 (Am 2 n - s)D A (Aml - s)D 5 (Aml - s)As(4m 2 - s) 



(38) 
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Li(s) 



(s - Am 2 n f + ai(s - Ami) 5 + a 2 (s - Ami) 4 + «s(s - 4m 2 ) 3 + 

+a 4 (s - Ami) 2 + a 5( s ~ ^ m l) + a 6+ 
+y/s(c 1 (s - Ami) 5 + c 2 (s - Ami) 4 + Cs(s - 4m 2 ) 3 + 

+ c 4 (s - 4m 2 ) 2 + c 5 (s - 4m 2 ) + c 6 



(39) 



A(s) = m 2 - s - ^11^ (s) , 



(40) 



16n 



ilRab 



P^{s) = P nl (s - ie) = P nl {s) + 2ip W7r {s) (l + 04^1 (s)) 



(41) 



(42) 



here 



A 2 _i_ e _ A-m 2 
4m£ 



(43) 



= 1 / 2 ^ 



1 r 00 s-4m 2 p^(s') 



4m 2 s' — s — 2£ 



D 1A (Am 2 w - s)D 2A (Aml - s) J D 3A (4m 2 - s) 



(44) 
(45) 



L 2 ( S )=4m 2 ( S 2 + /3 S + 7 i S 3/2 + 7 2 s 1/2 



(46) 



P* 2 {s) = P^{s - it) = P n2 (s) - 2ip n7T (s)k 2 K 2 (s) . 



(47) 



Note that this parametrization was inspired by Ref. 19[, devoted to proof that the 
propagators (jHJ) satisfy the Kallen-Lehmann representation in the wide domain of coupling 
constants of the scalar mesons with the two-particle states. Following the ideas of this paper 
the conditions 



K x (s) > 0, K 2 {s) > at s > Ami 



10 



guarantee absence of poles on the physical sheet in Eq. (135|) (of course, the restrictions 
of Sec. |IV] should be fulfilled too). Note also that we choose the denominator of fl35|) as 
P ff i(s)P 7r2 (s) for our comfort. 



IV. RESTRICTIONS ON THE PARAMETERS 



Some parameters are fixed by the requirement of the proper analytical continuation of 

q and s = m\ A 



amplitudes. The denominators P n i and P n 2 have zeroes at s = ml and s = m\ A respectively, 



both belonging to the interval < s < 4m^. These zeroes should be compensated by zeroes 
in any pair from P* Xl P* 2 and S® res . We choose 

S° res (ml A ) = 0, (48) 



see Eq. (EHD [28| 



The requirement of the Tq finiteness at s = leads to 2 conditions. Really, on the real 
axis for s > 4m^ we have 

Oback _ P:i(s) _ Pnl(s-ie) ^ l+JfrQQ 

1 ~ P^s)- P nl (s + ie)- i + ZiP ™ [S) P^s) ' 
.back _ p ^(s) _ PA S ~ ie) _ 1 .-_ ^(s) 



S-f cfc = = ^ = 1 - 2ip^(s) 



PAS) Pn2(s + ie) ^ KJ P,2(S)' 

So, to avoid singularity in the 



T, 



o 



Cjback gback gO res 



2ip^(s) 
at s = 0, where p^s) becomes infinite, we require 



l + ^i(0) = 0, 

as for ^(O), it is equal to zero at s = via construction, see Eq. (J4"5]) . Note that, 
alternatively, one may require Tq res (0) = 0. 



Additionally, as it was noted in Refs. 
equations imposes the condition 



29|], crossing symmetry implemented by Roy 
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^l( m 2 = 0) = 0. 
dm 

Recall that the condition Eq. ( l30l) removes the singularity in the T(tttt — > KK) ampli- 
tude. One can see that no special prerequisite to Adler zero existence in the nir scattering 
amplitude should be imposed, because it appears when we take into account the results of 
Ref. fl. 



V. DATA ANALYSIS 



Analyzing data, we imply a scenario motivated by the four-quark model 30J, that is, 
the cr(600) coupling with the KK channel, g a K+K-, is suppressed relatively to the coupling 
with the 7T7T channel, (y( (T7r + 7r -, the mass of the a meson m a is in the 500-700 MeV range. In 
addition, we have in mind the Adler self-consistency conditions for the Tq(tttt — > tttt) near 
;he 7T7T threshold. The general aim of this section is to demonstrate that the data and the 
6] results on the 7T7t amplitude are in excellent agreement with this general scenario. 

As in Ref. Q for — > 7r°7r°7 decay we use the KLOE data Si for m > 660 MeV. For 



the 5q we use the "old data" 



31 



35j, 44 points up to 1.2 GeV [36]. Besides, we take into 



37 



account the new precise data in the low energy region 

The inelasticity 77$ (m) and the phase 5 7TK (m) of the amplitude T{im — > KK) are essential 
in the fit region, 2m^+ < m < 1.2 GeV. As for the inelasticity, the experimental data of 



Ref. 



31] gives evidence in favor of low values of ^(m) near the KK threshold. At present 



the contribution of the 7777, 7/77, and 7/7/ channels does not affect much the overall fits. To 
fix a relation between the KK and 7777 channels reliably the inelasticity should be measured 
with accuracy 10 times better than the existing one. The situation with the experimental 
data on 8* K {m) is controversial and experiments have large errors. We consider these data 
as a guide, whose main role is to fix the sign between signal (jlj) and background amplitudes 



3]), and hold two points of the experiment 39], see Fig. |9j As for inelasticity, for fitting we 



used only the key experimental point rj^im = 1.01 GeV) = 0.41 ± 0.14, see Fig. |5j 

Providing all the above conditions, we have obtained perfect agreement with the general 
scenario under consideration, see Fits 1, 2 in Tables I, II, and III and Figs. 1-10. Fits 1 and 
2 show that the allowed range of cr (600) and /o(980) parameters is rather wide. 

The values of gj oK+K -/4:ii in Fits 1 and 2 (1 GeV 2 and 2 GeV 2 , correspondingly) show a 
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scale of possible deviation of this constant. This may be important to coordinate g"j oK+K - / 4tt 
with gl oK+K -/4:7i |40[], note the latter is usually larger than 1 GeV 2 . 

In addition, we carry out Fit 3, where cr(600) and / (980) are coupled only with the Tin 
channel. As seen from Table I and Figs. 11-13, Fit 3 is in excellent agreement with the data 
on the Sq up to 1 GeV and the [6J results. 

We introduce 52 parameters, but for restrictions (expresses 5 parameters through others) 
and parameters (or their combinations), that go to the bound of the permitted range (7 
effective links), the effective number of free parameters is reduced to 40. It is significant 
that fits describe not only the experimental data (about 80 points), but also the mr amplitude 
from the 6| in the range —5m 2 < s < 0.64 GeV 2 which is treated along with experimental 
data. 



, see Table I. One 
6| approach. The 



The cx(600) pole positions, obtained in Fits 1 and 2, lie far from Eq. ([T] 
of the possible reasons is neglecting KK and other high channels in the 
role of high channels can be estimated with the help of Fit 3, whose <r(600) pole position is 
considerably closer to Eq. (CQ), see Table I. 

Note that kernels of the background integrals ( 1381) and ( ]4"5i) are positive in the range of 
integration [2m„.,oo), Fit 1 kernels are presented in Fig. 7. 

The Adler zero in the Tq (tttt — > tttt) is near s = (100 MeV) 2 in all Fits because we describe 
the amplitude fl. Fit 2 also has Adler zero in the T(tttt KK) at s = (166 MeV) 2 , Fit 1 
has a zero in the T(tttt KK) at s = -(601 MeV) 2 . 

The resonance amplitude Tq res has poles on the unphysical sheets of its Riemannian 
surface. As we have a multichannel case, the amplitude has the set of lists depending on 
lists of the polarization operators IT^(s). We show resonance poles only on some lists, see 
Tables IV and V. For this choice, in case of metastable states, decaying to several channels, 
the imaginary parts of pole positions M R would be connected to the full widths of the 
resonances [2ImM# = Yr = J2abF(R —> ab)]. Note that cr(600) and / (980) poles, shown in 
Table I, correspond to the first lines of Tables IV and V. 

As to the background amplitude T^ back , it has poles on the second sheet of the Rieman- 
nian surface, where P n i = or P n2 = 0. The P n i has a zero at s = (1246 — 104 i) 2 MeV 2 
for Fit 1, at s = (1354 - 110 if MeV 2 for Fit 2, and at s = (1056 - 142 i) 2 MeV 2 for Fit 3. 
The P n2 has a zero at s = (0.2 - 9.5 i) m 2 for Fit 1, at s = (2.0 - 8.9 z) m 2 for Fit 2, and at 
s = (—0.6 — 8.6 i)m 2 for Fit 3. These poles lie outside of the region 4m 2 < Re(s) < (1.2 
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GeV) 2 except the pole at s = (1056 - 142 if MeV 2 for Fit 3, but for this Fit the upper 
bound is 1 GeV 2 . 



Table I. Properties of the resonances and main characteristics are shown. 



r it 


1 
1 


o 
Z 


o 
O 


„ TV J \ T 

rrif , MeV 


979.1b 


nop crn 


9o4.(Jl 


9f K+K~, ^ev 


3.54 


5.01 





a' 2 

f0 4n , GeV 2 


1 


2 





9f n+n-, GeV 


-1.3737 


-2.1185 


0.3183 


Jbg^ GeV 2 


0.150 


0.357 


0.008 


x fo 


0.6640 


0.9584 


- 


r /o (m /o ),MeV 


55.2 


130.3 


3.0 


/o(980) pole, MeV 


986.2 - 25.5 i 


990.5 - 19.4 i 


978.9 - UAi 


m ff , MeV 


487.59 


506.95 


480.46 


SW+tt-, GeV 


2.7368 


2.6735 


2.5871 


GeV 2 


0.596 


0.569 


0.533 


9aK+K-, GeV 


0.552 


0.774 





GeV 2 


0.024 


0.048 







0.9750 


0.8201 





r.K), MeV 


377.8 


352.9 


340.2 


<j(600) pole, MeV 


581.0 - 212.7 i 


613.8- 221.4 % 


528.6 - 220.3 i 


C, GeV 2 


0.04317 


-0.07633 


-0.11734 




-70.62 


-73.6 


- 


mi, MeV 


801.90 


814.88 


- 


m 2 , MeV 


465.95 


554.95 




A x , GeV 


1.142 


1.030 




a° , m- 1 


0.223 


0.226 


0.221 


Adler zero in tttt — > nn 


(94.4 MeV) 2 


(96.8 MeV) 2 


(87.1 MeV) 2 


77^(1010 MeV) 


0.55 


0.45 




X 2 ph ase ( 44 Points) 


45.9 


50.6 


26.3 (34 points) 


X% (18 points) 


24.9 


19.1 
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Table II. Parameters of the first background (P n i) are shown. 



Fit 


1 


2 


3 




-3.105 


-4.549 


-1.498 


a 2 


0.01136 


0.00998 


0.05821 


a 3 











a 4 


4.9328 


13.1111 


1.2475 


an, GeV 2 


604.137 


624.512 


-792.804 


a 2 , GeV 4 


920.111 


1000.739 


-384.477 


a 3 , GeV 6 


785.958 


781.770 


416.645 


a 4 , GeV 8 


223.623 


211.195 


198.772 


a 5 , GeV 10 


24.5339 


23.8517 


25.4265 


a 6 , GeV 12 


0.248657 


0.314094 


0.198560 


ci, GeV 


356.128 


224.404 


995.905 


c 2 , GeV 3 


-2735.40 


-2600.82 


-1070.75 


c 3 , GeV 5 


284.008 


445.192 


542.745 


c 4 , GeV 7 


430.758 


461.717 


411.927 


c 5 , GeV 9 


49.7913 


47.2357 


51.4206 


c 6 , GeV 11 


-0.664290 


-0.684002 


-0.635647 


mi, MeV 


1105.67 


1111.87 


1002.31 


9i, MeV 


347.70 


350.48 


306.18 


m 2 , MeV 


1061.53 


1141.92 


806.93 


92, MeV 


344.12 


381.73 


350.51 


m 3 , MeV 


1061.85 


1169.51 


781.76 


<?3, MeV 


311.56 


311.80 


322.57 


m 4 , MeV 


970.78 


1040.96 


970.78 


<?4, MeV 


457.52 


455.56 


376.88 


m 5 , MeV 


1176.39 


1320.55 


1153.21 


05, MeV 


544.43 


588.48 


500.59 


m 6 , MeV 


1521.20 


1621.10 


1808.74 


9e, MeV 


739.93 


750.75 


841.57 
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Table III. Parameters of the second background (P n 2) are shown. 



Fit 


1 


2 


3 


A, MeV 


83.238 


74.477 


70.268 


k 2 


0.0152934 


0.0168176 


0.0150655 




239.184 


221.055 


263.511 


7i 


1006.367 


928.743 


878.056 


72 


22.7004 


23.3341 


29.4097 


m 1A , MeV 


491.92 


84.77 


687.43 


9ia, MeV 


469.29 


492.03 


364.68 


m 2 A, MeV 


531.81 


639.95 


528.40 


92A, MeV 


452.20 


261.48 


378.65 


m 3 A, MeV 


670.64 


565.16 


608.72 


93A, MeV 


299.23 


428.97 


370.98 



Table IV. o~(600) poles (MeV) on different sheets of the complex s 
plane depending on lists of polarization operators U ab (s) are shown 



U KR list 


IT" list 


n^' list 


IT^' list 


Fit 1 


Fit 2 


I 


I 


i 


I 


581.0- 212.7 i 


613.8 - 221.4i 


II 


I 


i 


I 


617.5 -353.0 i 


609.8 — 291.6z 


II 


II 


i 


I 


554.3 -375.3 i 


559.4- 346.6 i 


II 


II 


ii 


I 


579.0 - 475.2 i 


569.7 - 410.7 i 


II 


II 


ii 


II 


625.7- 474.9 i 


581.6 - 411.0i 



Table V. / (980) poles (MeV) on different sheets of the complex s 
plane depending on lists of polarization operators U ab (s) are shown 



U KR list 


n w list 


IT™' list 


IF'"' list 


Fit 1 


Fit 2 


I 


I 


I 




986.2 - 25.5 i 


990.5 - 19.4 i 


II 


I 


I 




916.9 - 299.4i 


1183.2 - 518.6i 


II 


II 


I 




966.8 - 450.5 i 


1366.0 - 756.5 i 


II 


II 


II 




962.6 - 465.2 i 


1390.7- 813.02 


II 


II 


II 


II 


962.5 - 608.0 % 


1495.6 - 1057.72 
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FIG. 1: The 7r°7r° spectrum in the <p 7r°7r°7 decay, theoretical curve, and the KLOE data (points) 
are shown: a) Fit 1, b) Fit 2. 




0-4 0.6 0.8 1 1.2 0.4 0.6 0.8 1 1.2 

m, GeV m# GeV 



(a) (b) 

FIG. 2: The phase 5q of the nir scattering (degrees) is shown: a) Fit 1, b) Fit 2. 

Note it would be naive to treat the poles in the background as resonances [/o(1370), for 
example] because in our approach to consider additional resonances one should extend the 
matrix of the inverse propagators, etc. 



VI. CONCLUSION 



Thus, the background phase ( 1341) allows us to obtain proper analytical features of the 
7T7T scattering amplitude, link results of js] with properties of light scalars simultaneously 
describing experimental data. The obtained description is in agreement with the scenario 
based on the four-quark model. The main signatures of this scenario are the weak coupling 
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280 300 320 340 360 380 400 
m, MeV 



(a) 



280 300 320 340 360 380 400 
m, MeV 



(b) 



FIG. 3: The comparison of the experimental data on 5® [3j| and the obtained curve is shown: a) 
Fit 1, b) Fit 2. 
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280 300 320 340 360 380 400 420 
m, MeV 




280 300 320 340 360 380 400 420 
m, MeV 



(a) 



(b) 



FIG. 4: The comparison of the experimental data on <5[j [381 ] and the obtained curve is shown: a) 
Fit 1, b) Fit 2. 

of the <t(600) meson with the KK channel compared to the ttti one and the weak coupling 
of the /o(980) meson with the tttt channel compared to the KK one, see Table I, that results 
in the weak cr(600) — /o(980) mixing |2j. The ratios (g a K+K- / 9an+n-) 2 ~ 0.05 — 0.1 and 
{gfo-rr+ir- 1 9f Q K+K-) 2 ~ 0.15 — 0.18, see Table I, indicate roughly that the 90 percentage of 
cr(600) is udud and the 80 percentage of /o(980) is sdds. 

Resonance masses and widths thr and r^(m^) in our formulas (which may be called 
" Breit- Wigner" masses and widths) have clear physical meaning, in contrast to the resonance 
poles in the complex plane. At first, what sheet of the complex plane should be considered? 
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m, GeV 
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(b) 
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FIG. 5: The inelasticity t]q is shown: a) Fit 1, b) Fit 2. 




500 600 
m, MeV 



(a) 



700 



800 



80- 
60 
40 
20 




(b) 



300 400 500 600 700 800 
m, MeV 



FIG. 6: The phase <5g of the 7T7T scattering is shown. The solid line is our description, dashed lines 
mark borders of the corridor [6||, and points are experimental data: a) Fit 1, b) Fit 2. 

For a (600) it is natural to consider the first line of Table IV [at any rate, it would be 
correct for very narrow cr (600)]. The obtained pole positions in this case do not agree with 
the pole position obtained in Ref. jfj], see Eq. ([I]). Note that the a(600) pole position 
is dictated by the <r(600) propagator in our case, because the <r(600) — /o(980) mixing 
is small. Providing the pole position ([1]) and taking into account only the tttt channel 
in the propagator, we can determine <r(600) mass and coupling to the tttt channel, and the 
obtained values contradict the Kallen-Lehmann representation, see 19] . Taking into account 
additional channels we may fulfill the Kallen-Lehmann representation, but the region of 
permitted <r(600) parameters do not allow us to describe experimental data in the current 
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-4 



(a) 



(b) 



FIG. 7: The real and the imaginary parts of the amplitude Tq of the mr scattering (s in units 
of m 2 ) are shown. Solid lines show our description, dashed lines mark borders of the real part 
corridor and the imaginary part for s < [6|: a) Fit 1; b) Fit 2. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 

m, GeV 



(a) 



50 








40 






30 






20 






10 




A 





20 30 
m, GeV 



(b) 




(c) 



FIG. 8: Kernels K^m 2 ) and K 2 (m 2 ) for Fit 1 are shown: a) Ki(m 2 ) below 2 GeV. The minimum 
near 1.4 GeV is 0.25. b) Ki(m 2 ) up to 50 GeV, then it asymptotically tends to 1. c) K2(m 2 ) up 
to 2 GeV, then it asymptotically tends to zero. 



model. 

Note that the Roy equations are approximate, they take into account only the irn channel. 
This can lead to a different analytical continuation and, hence, explain deviation of the 
cr(600) pole position, compare Fit 3 with Fits 1 and 2 in Table I 41]. 

The current activity, aiming extremely precise determination of the cx(600) pole position, 
has taken the forms of the Swift's grotesque. Really, the residue of the a pole can not be 
connected to coupling constant in the Hermitian (or quasi- Hermitian) Hamiltonian, see Ref. 
[5I, for it has a large imaginary part and this pole can not be interpreted as a physical state 
for its huge width. 
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FIG. 10: The \P K (m)\ 2 is shown, see Eq. (7): a) Fit 1; b) Fit 2. 

The futility of the approach that is based on the poles treatment may be additionally 
illustrated by Fit 2. As seen on line 1 of Table V, the real part of the / (980) pole ReMf on 
the II sheet of the T ° exceeds the K + K~ threshold (987.4 MeV), it means that ImM fo equals 
to -(r(/ (980) ->■ tttt) - r(/ (980) ->■ K + K')^j /2, which is physically meaningless. In this 
case we should take Yl K+K from the second sheet, this gives the pole at Mf = 989.6 — 168.7 i 
MeV, with ReM fo between the K + K~ and K°K° thresholds again. As we work on the s 
plane, we should consider not Mf , but Mj Q = (0.951 — 0.334 i) GeV 2 . So, we have the pole 
with a real part below the K + K~ and K°K° thresholds and an imaginary part dictated by 
analytical continuation of the kaon polarization operators. 

To reduce an effect of heavier isosinglet scalars we restrict ourselves to the analysis of the 
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FIG. 11: The phase <5g of the irir scattering, Fit 3 and the experimental data are shown. 
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280 300 320 340 360 380 400 420 
m, MeV 



(b) 



FIG. 12: The phase <5§ of the irir scattering, Fit 3 is shown. The comparison with the data is 



available from a) [37[, b) 



n 



mass region m < 1.2 GeV. As to mixing light and heavier isosinglet scalars, this question 
could not be resolved once and for all at present, in particular, because their properties are 
not well established up to now. A preliminary consideration was carried out in Ref. 42), 
where, in particular, it was shown that the mixing could affect the mass difference of the 
isoscalar and isovector. 

The factor |P^(s)| 2 modifying the \T(ttti — > KK)\ 2 , see Eqs. (7) and (l3Tj) . is shown on 
Fig. [TOJ This factor does not change the kaon loop model radically, but helps to fulfill the 
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FIG. 13: a) The phase 5% of the 7T7T scattering, Fit 3 is shown. The solid line is our description, 
dashed lines mark borders of the corridor [6j|, and points are experimental data, b) The real and the 
imaginary parts of the amplitude Tq of the mr scattering (s in units of m^) are shown. The solid 
lines correspond to Fit 3, dashed lines mark borders of the real part corridor and the imaginary 
part for s < [g]. 

requirement f[3"Uj) and to improve the data description. The influence of this factor may be 
reduced in order to use a more skillful form than Eq. (l3Tj) for it. 

New precise experimental data are needed for the investigation of light scalars. The 
elucidation of the situation, a contraction of the possible variants or even the selection of 
the unique variant, requires considerable effort. The new precise experiment on Tin — > KK 
would give the crucial information about the inelasticity 77$ and about the phase 8g K (m) 
near the KK threshold. The forthcoming precise experiment in KLOE on the — > 7r°7r°7 
decay will also help to judge this phase in an indirect way. The precise measurement of the 
inelasticity rj® near 1 GeV in tttt — > Tin would also be very important. 



It is of interest to update our analysis of the —> a Q (980)7 ~^ V^ ! decay 23] and 
the 77 — > ao (980) — > -k°7] description |43] in this analytical approach. Probably, such an 
approach would also be useful for the k(900) meson investigation in the irK channel. 
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